arXiv:nucl-th/0112083vl 31 Dec 2001 


Kaon and Pion Fluctuations from Small Disoriented Chiral Condensates 


Sean Gavin“ and Joseph 1. Kapusta*' 

a) Department of Physics and Astronomy, Wayne State University, Detroit, MI, 48202 

b) School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455 

(February 8, 2008) 

Enhancement of D, and baryon production in Pb+Pb collisions at a c.m. energy of 17 A GeV can be 
explained by the formation of many small disoriented chiral condensate regions. This explanation 
implies that neutral and charged kaons as well as pions must exhibit novel isospin fluctuations. 
We compute the distribution of the fraction of neutral pions and kaons from such regions. We 
then propose robust statistical observables that can be used to extract the novel fluctuations from 
background contributions in and KgK^ measurements at RHIC and LHC. 


I. INTRODUCTION 

Heavy ion collisions at the Brookhaven Relativistic 
Heavy Ion Collider (RHIC) at center of mass energies 
up to 200 A CeV and the CERN Large Hadron Collider 
(LHC) at 5.5 A TeV may produce matter in which chiral 
symmetry is restored. One possible consequence of the 
restoration and the subsequent re-breaking of chiral sym¬ 
metry is the formation of disoriented chiral condensates 
(DCC) - transient regions in which the average chiral or¬ 
der parameter differs from its value in the surrounding 
vacuum tB _ 

Measurements of ft and Ut baryon enhancement at 
17 A CeV at the CERN SPS can be explained by the 
production of many small DCC regions within individual 
collision events H] . If true, this explanation has two im¬ 
portant consequences. First, the DCC regions must be 
rather small, with a size of about 2 fm. Such a size is con¬ 
sistent with predictions based on dynamical simulations 
of the two flavor linear sigma model |^. More startling 
is the second implication that the evolution of the con¬ 
densate can have a significant effect on strange particle 
production. The importance of strange degrees of free¬ 
dom in describing chiral restoration has been long appre¬ 
ciated |)7[-pl|, but simulations of the three flavor linear 
sigma model had suggested that strange kaon fields are 
much less important than the pion fields |^. Neverthe¬ 
less, the Ul and H data demand that we explore without 
prejudice techniques for measuring kaon fluctuations. 

In this paper we study pion and kaon isospin fluctua¬ 
tions in the presence of many small DCC. In the next sec¬ 
tion we compute probability distributions that describe 
the DCC contribution to these fluctuations. Pion fluctu¬ 
ations due to many small DCC have been addressed by 
Amado and Lu [|^ and Chow and Cohen Q, although 
the distribution we compute is new. Ours is the first 
work to study kaon fluctuations. In sec. 3 we combine 
the DCC fluctuations with a contribution from a random 
thermal background. In sec. 4 we discuss how the size 
and number of DCC vary with impact parameter, target 
and projectile size. In sec. 5 we assess robust statistical 
observables that can be used to measure the impact of 


many small DCC at RHIC and LHC. In particular, we 
obtain a dynamic isospin fluctuation observable analo¬ 
gous to the dynamic charge observable used to measure 
net charge fluctuations at RHIC Of the quantities 
considered, this observable isolates the DCC effect from 
other sources of fluctuations best. 

To illustrate how a strange DCC can form, we first con¬ 
sider QCD with only up and down quark flavors. Equi¬ 
librium high temperature QCD respects chiral symmetry 
if the quarks are taken to be massless. This symmetry 
is broken below Tc ~ 150 MeV by the formation of a 
chiral condensate (cr) ~ Ifau + dd) that is a scalar isopin 
singlet. However, chiral symmetry implies that a is de¬ 
generate with a pseudoscalar isospin triplet of fields with 
the same quantum numbers as the pions. In reality, chi¬ 
ral symmetry is only approximate and the 140 MeV pion 
mass is different from the 800 ± 400 MeV mass of the 
leading sigma candidate [|^ . Nevertheless, lattice calcu¬ 
lations exhibit a dramatic drop of (cr) near Tc at finite 
quark masses. 

A DCC can form when a heavy ion collision produces a 
high energy density quark-gluon system that then rapidly 
expands and cools through the critical temperature. Such 
a system can initially break chiral symmetry along one 
of the pion directions, but must then evolve to the T = 
0 vacuum by radiating pions. A single coherent DCC 
radiates a fraction of neutral pions compared to the 
total that satisfies the probability distribution 

= 0<U<1, ( 1 ) 

||r^-[l^ . Such isospin fluctuations constitute the primary 
signal for DCC formation. The enhancement of baryon- 
antibaryon pair production is a secondary effect due to 
the relation between baryon number and the topology of 
the pion condensate field [p0| . 

This two flavor idealization only applies if the strange 
quark mass rus can be taken to be infinite. Alternatively, 
if we take rUg = = rud = 0, then the chiral conden¬ 

sate would be an up-down-strange symmetric scalar field. 
The more realistic case of TOs ~ 100 MeV is between these 
extremes, so that (cr) ^ {cos 9{uu + dd) -I-sin0(ss)). The 
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mixing angle 9 is highly uncertain since it depends on 
the sigma mass together with the tt^ K,r] and t]' masses 
and the r] — rj' mixing angle . A disoriented conden¬ 
sate can evolve by radiating tt, K, rj and 77 ' mesons, with 
the neutral pion fraction satisfying (1). Randrup and 
Schaffner-Bielich find that the kaon fluctuations from a 
single large DCC satisfy 


Limit Theorem. In the next section these distributions 
will be folded together with a random or thermal source 
which most likely would comprise the bulk of the mesons 
in a high energy heavy ion collision. 


A. Kaons 


Piifx) = 1 0 < /x < 1, (2) 

where Jk = {K^+T^)/{K^ +K~ +K^+K^). Moreover, 
the condensate fluctuations can now produce strange 
baryon pairs |^. Linear sigma model simulations indi¬ 
cate that pion fluctuations dominate three-flavor DCC 
behavior, while the fraction of energy imparted to kaon 
fluctuations is very small due to the kaons’ larger mass. 
On the other hand, domain formation may be induced by 
other mechanisms such as bubble formation or decay 
of the Polyakov loop condensate . 

Why does the DCC’s size matter? Pion measurements 
in individual collision events can distinguish DCC isospin 
fluctuations from a thermal background only if the dis¬ 
oriented region is sufficiently large |^. DCC can then 
be the dominant source of pions at low transverse mo¬ 
menta, since (pi) ~ 1/i? for a coherent region of size R. 
Experiments focusing on low pt can study neutral and 
charged pion fluctuations [|^ , wavelet and HBT sig¬ 
nals to extract detailed information. In contrast, 

for small domains {R < 3 fm |^) DCC signals are hid¬ 
den by fluctuations due to ordinary incoherent produc¬ 
tion mechanisms. This holds even if many such regions 
are produced per event. DCC mesons from small regions 
may have momenta of a few hundred MeV, nearer the pp 
mean value. Different regions would not add coherently 
to alter HBT, nor would their small spatial structures 
affect wavelet analyses. 

Importantly, baryon pair enhancement is substantial 
only if there are many small incoherent regions. The large 
winding numbers that produce baryon-antibaryon pairs 
require many small regions with random relative orienta¬ 
tions of the pion held. To describe strange antibaryon en¬ 
hancement, Kapusta and Wong assume roughly 100 DCC 
regions of size of roughly 2 fm |^ . Topological models of 
baryon-antibaryon pair production successfully describe 
e -I- e— and hadronic collision data . The connection 
of DCC to topological pair production was pointed out 
in Ref. ||^; see also p^ . 


II. FLUCTUATIONS IN NEUTRAL DCC 
MESONS 


Dehne f = {K° + K°)/{K+ + R- + K° + Tf). To 
an excellent approximation the number of neutral kaons 
is equal to twice the number of short-lived neutral kaons 
Ks which are more readily measurable in high energy 
heavy ion collisions. The fraction / ranges from 0 to 1. 

The statistical distribution in / for a single domain is 
pi(/) = 1. The distribution for n randomly oriented, 
independent domains is 


Pnif) = / n d/fc Pl(/fc) W / - - E I ■ (3) 

The Dirac delta function can be represented as an inte¬ 
gral. Then can be written as 


Pn{f) = 


dz 


-ifz 


^ixz / n 


dxpi{x)e^ 


( 4 ) 


Since pi{f) = 1, the integration over x can be done, 
resulting in a one- dimensional integral. 

= cos [n(l - 2/)t] (5) 

The integral can be evaluated and expressed in terms of 
a hnite sum. 

(>,.(/) ="’ E (-ij'hAALJTl (6) 


0<fe<n(l —/) 


fc!(n — fc)! 


It is useful to have a simple analytic formula for pn in 
the limit that n ^ 1. In this limit the factor (sint/t)" 
in the integral formula is strongly peaked at t = 0. Let 
us write this factor as exp(—E(t)) with a view towards a 
saddle point approximation. We get 

F(0) = F'(0) = 0 
F"(0) = n/3 

. (7) 


Use of this approximation yields the asymptotic formula 


In this section we will compute the statistical distri¬ 
bution of the ratio of neutral to total number of mesons, 
first for kaons and then for pions. In both cases the limit 
that the number of DCC domains becomes large is taken. 
It is natural that this limit results in a Gaussian distribu¬ 
tion for both kaons and pions on account of the Central 


Pn(/) = \/—exp[- 6 n(/- 1 / 2 ) 2 ]. 

V TT 

The distribution is strongly peaked around / = 1/2 as 
one might expect. 
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FIG. 1. The probability distribution for the ratio of neutral 
to total number of DCC kaons from n domains. The dashed 
curves represent Gaussian distributions. From eqs. (6) and 
( 8 ). 

Figure 1 shows the evolution of Pn{f ) with n. It goes 
from a flat distribution for n = 1 to a Gaussian sharply 
peaked at / = 1/2 as n becomes large compared to 1. In 
fact a Gaussian is a very good representation for n > 2. 


B. Pions 


Define / = 7 r°/( 7 r+ + 7 r“ + 7 r°). To a good approxi¬ 
mation the number of neutral pions is equal to half the 
number of photons. Therefore, to this level of precision, 
it is not necessary to identify each 7 r° via its decay into 
27 . The fraction / ranges from 0 to 1. 

The statistical distribution in / for a single domain is 
Pi{f) = The distribution for n randomly ori¬ 

ented, independent domains can be computed along the 
same lines as for kaons. 


Pnif) 


r dz 

f dx ixzjn 

/-oo 27r 

Jo 2y/x 


(9) 


Since pions from DGG have been extensively studied we 
shall be content to evaluate the distribution in the large n 
limit. This is accomplished by expanding the exponential 
in the x integration to second order in 1 /n, evaluating the 
resulting integrals, and exponentiating. Thus 


dx 

2^/x" 


^ixz / n 


exp 


3n 


2^2 1 
45n2 


( 10 ) 


The z integral can then be done, yielding a Gaussian 
centered at / = 1/3. 


Pnif) 



exp —45n(/ — 


1/3)V4] 


( 11 ) 


III. FOLDING DCC AND THERMAL MESONS 


In a more realistic scenario some kaons will come from 
the decay or realignment of DCC domains and some will 
come from more conventional sources. We shall refer to 
the latter as random or thermal, even though that may 
be a bit of a misnomer. What we mean by random or 
thermal is that the distribution of kaons from non-DCC 
sources is 

Poifo) = exp [-(/o - Il2)^l2al] . (12) 

For a completely random source the width ctq is related 
to the total number Arandom of non-DCC kaons by 

^ 1/2(1-1/2) ^ 1 ^ 

Arandom 4A,.andom 


Now let US assume that a fraction uk of all kaons come 
from non-DCC sources and the remaining fraction Pk = 
1 — ax come from n ^ 1 independent DCC domains. 
Letting A denote the total number of kaons, we have 
Arandom = Oic A and Adcc = /3ic A. Folding together 
two Gaussians gives a Gaussian. 


PkU) = J dfodfnPoifo)Pnifn)Sif - axfo - 13k fn) 

exp[-(/-l/2)V2Ay (14) 


\/27rA|. 

The net width is 


A2 - 


OtK 

4A 


12n 


4A 


+ 


/3|: _ ^ 

12n 4 A 


(15) 


The expression in curly brackets at the end represents 
the difference between the actual width and the width the 
distribution would have if there was no contribution from 
DCC kaons. This change in the width may be positive 
or negative, depending on the parameters. 

An analogous analysis can be given for pions. This 
results in the distribution 


pM) 


J dfodf„poifo)Pnifn)Sif - a^fo - Pnfn) 

exp [-(/ - 1/3)V2A2] (16) 


with a net width of 


9A 45n 9A ( 45n 9A j 


(17) 


As with the kaons, the last expression in curly brackets 
represents the difference between the actual width and 
the width the distribution would have if there was no 
contribution from DCC pions. Note that the fractions 
ax and need not be the same. 
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IV. VOLUME OR SURFACE SCALING? 

The issue we wish to address is whether the number 
of DCC mesons (kaons or pions) scales with the volume 
or surface area of the system. This is an important issue 
when studying the impact parameter dependence or the 
dependence on the size of the projectile and target nuclei. 

In this paper it is assumed that DCC have a typical 
size of order 2 fm which does not change much with col¬ 
lision energy or the total volume of the system. Thus 
the number of domains n is just given by the ratio of the 
two volumes: n = I4ystem/TDCC- Scaling of the num¬ 
ber of DCC kaons or pions, Nucc, with n or I4ystem is 
the same because the size of individual domains is fixed. 
The Njjcc will depend on the extra energy associated 
with the formation of a domain. If the up and down 
quark masses were zero then QCD would have perfect 
SU(2) flavor symmetry. In that case the energy density 
of a large uniform domain would be independent of its 
orientation. All directions in chiral space are equivalent. 
However, the misalignment between adjacent domains re¬ 
sults in a surface energy, and so the number of DCC pions 
would be proportional to the total surface energy between 
domains. The up and down quark masses are not zero 
(they’re about 5 to 7 MeV), and this will result in an ex¬ 
cess volume energy too. For pions one might expect the 
surface energy to dominate since the up and down quark 
masses are so small. For kaons one might expect the vol¬ 
ume and surface energy contributions to be comparable 
on account of relatively large mass of the strange quark 
(about 120 MeV). 

Let us analyze the scaling of the excess surface en¬ 
ergy more quantitatively. Consider a cube with sides of 
length L into which ht n = {L/Vf cubic domains, each 
with sides of length 1. Assuming that each domain is 
oriented independently of its neighbors, the total sur¬ 
face energy scales with the total surface area, which is 
3 (j -f l) With L ^ I, and with the domain size I 
fixed, the total surface area, energy, and therefore num¬ 
ber of DCC mesons scale with n to the power of one. 
Thus Nucc oc n (X I4ystem no matter whether one imag¬ 
ines the excess energy being associated with domain in¬ 
terfaces or with domain interiors. 


V. STATISTICAL ANALYSIS 

Detection of small incoherent DCC regions in high en¬ 
ergy heavy ion collisions requires a statistical analysis in 
the TT^TT^ or the KgK^ channels. Neutral mesons can 
be detected by the decays 7 r° ^ 77 or Kg tt+tt”. 
The analysis we propose is sensitive to correlations due 
to isospin fluctuations. We expect these correlations to 
vary when DCC regions increase in abundance or size 
as centrality, ion-mass number A, or beam energy are 
changed. Correlation results combined with other sig¬ 
nals, such as baryon enhancement [^, can be used to 


build a circumstantial case for DCC production. 

Correlations of 7 r° 7 r=*= and KgK^ can be determined 
by measuring the robust isospin covariance, 


Rco 


(iVeiVo) - (W)(jVo) 
{Nc){No) 


(18) 


where Nq and are the number of neutral and charged 
mesons. We take Nq = N^o and Nf. = N^++N^- for pion 
fluctuations and Nq = 2Nj^o and Nc = Nx+ + N^- for 
kaon fluctuations. The ratio (18) has two features that 
are convenient for experimental determination. First, 
this observable is independent of detection efficiency as 
are the “robust” ratios discussed in . Robust observ¬ 
ables are useful for DCC studies because charged and 
neutral particles are identified using very different tech¬ 
niques and, consequently, are detected with different ef¬ 
ficiency. Observe that robust quantities are not affected 
by the unobserved AT®, since the strong-interaction eigen¬ 
states and are a superposition and K'g un¬ 
til their decay well outside the collision region. Second, 
since (18) is obtained from a statistical analysis, individ¬ 
ual 7 r° ^ 77 or ATg ^ 7 r“'" 7 r“ need not be fully recon¬ 
structed in each event. This feature is crucial because it 
would be extraordinarily difficult - if not impossible - to 
reconstruct a low momentum 7 r° in heavy ion collisions 
except on a statistical basis. 

Next we define robust variance 


{Nl) - {Naf - {No) 
{Na? 


(19) 


where a = c or 0. To see why (19) is robust, denote the 
probability of detecting each meson e and the probability 
of missing it 1 —e. For a binomial distribution the average 
number of measured particles is {NaY^'^ = e(A^a) while 
the average square is {N^Y^^ = + ^(1 ~ Yi^a)- 

We then find 


= Raa, ( 20 ) 

independent of e the proof that (18) is robust is 
similar. The ratios (18) and (19) are strictly robust only 
if the efficiency e is independent of multiplicity. Further 
properties and advantages of these and similar quantities 
are discussed in ||^ . 

To study DCC fluctuations we define the dynamic 
isospin observable 

^dyn = Rcc + Roo ~ 2i?c0- (21) 

Analogous observables have been employed to study net 
charge fluctuations in particle physic s and were 

considered in a heavy ion context in |ll5|] and [^l| . This 
quantity can be written in terms of 


V 


cO 




( 22 ) 
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To isolate the dynamical isospin fluctuations from other 
sources of fluctuations, one obtains (21) by subtracting 
from (22) the uncorrelated Poisson limit = (-^o)~^ + 
{Nc)~^. Indeed, we show in (27) below that the quantity 
(21) depends primarily on the fluctuations of the neutral 
fraction /, while the individual ratios (18) and (19) have 
additional contributions. 

We illustrate the effect of DCC on the dynamic isospin 
fluctuations by writing Nq = fN and = {1 — f)N. 
Small fluctuations on f or N results in the changes 


ANo 

(No) 

AW 

m 


AN Af 

M ^ 

AN Af 

M" 1^' 


We obtain the average 

{AN§) 2c A^ 

(WP {N){f) ^ {fy 


(23) 


(24) 


How big is the DCC effect compared to alternative 
sources of fluctuations? In the absence of DCC a = 1 
and /3 = 0 so that (29) implies = 0 for both pi- 
ons and kaons. On the other hand, in models which 
treat nuclear collisions as a superposition of indepen¬ 
dent nucleon-nucleon collisions, each nucleon-nucleon col¬ 
lision contributes an amount to the overall fluctu¬ 
ations. Consequently, M nucleon-nucleon collisions can 
contribute an amount /M to the total Q . While 
little is known from pp experiments about kaon fluctua¬ 
tions, HIJING and RQMD models yield negative values 
of For kaons, HIJING simulations of central Au-I-Au 
at 200 A GeV yield « -0.002 for 47 K+ and 44 
on average . The onset of a DCC contribution to ^dyn 
can substantially change this value. A detailed of analy¬ 
sis of this problem within microscopic models will appear 
elsewhere |]3^ . 


VI. DISCUSSION AND CONCLUSION 


Here the contribution of the variance of the total num¬ 
ber of mesons is u = (AV^)/(V)^ and the charge-total 
covariance is c = {ANAf). DCC formation primarily ef¬ 
fects the charge fluctuation contribution, A^ = ((A/)^), 
from (15) or (17). Similarly, 


(ANf) 2c A2 

7i^-^"(iV)(l-(/)) + (l-(/))2’ 

and 


RcO = V + 



l-(/) 


c 


A2 

( 1 -(/))" 


where Wo is given by (18). Using (21) we get 


dyn 


1 / A2 

(/)(!-(/)) l(/)(l-(/)) 



(25) 

(26) 

(27) 


This observable isolates the isospin fluctuations, whereas 
the individual Rat depend on the fluctuations in total 
meson number, v and c as well. 

We estimate the effect of DCC on the dynamical fluc¬ 
tuations (27) using (15) and (17). We take (N) — Nk 
for kaons and {N) = AU for pions; these are the total 
number of mesons of the indicated kind. For kaons 


vfy^iK DCC) = 4W . (28) 

and for pions 

ufy^irr DCC) = 4.5/3. • (29) 


These quantities can be positive or negative depending on 
the magnitude of j3 compared to the number of domains 
per kaon. In fact the dynamical fluctuation may even be 
positive for one kind of meson and negative for the other. 


Reference 1^ argued that the anomalous abundance 
and transverse momentum distributions of D and H 
baryons in central collisions between Pb nuclei at 17 A 
GeV at the CERN SPS is evidence that they are pro¬ 
duced as topological defects arising from the formation 
of many domains of disoriented chiral condensates (DCC) 
with an average domain size of about 2 fm. Motivated 
by this interpretation, we have studied the effect of DCC 
on the distribution of the fractions of neutral kaons and 
pions. We showed that the distributions are accurately 
described by Gaussians with centroids at / = 1/2 and 
1/3, respectively, once the number of domains exceeds 
just a few. Folding together kaons or pions arising from 
DCC with other sources that are Gaussian distributed re¬ 
sults once again in Gaussians. These may have a width 
that is greater or less than a purely random source with¬ 
out DCC formation. 

The DCC pioneers had hoped that a large 

percentage of pions might be emitted from just a few big 
domains, on the order of 5 to 8 fm (kaons were not con¬ 
sidered). Such large domains have been ruled out at SPS 
[^, but remain possible at RHIC. More conservatively, 
as the number of domains grow and their average size 
diminishes, the impression left on the fluctuations in the 
neutral fraction becomes more subtle and less unique. 
For many small domains, statistical measurements of 
both neutral kaons (pions) and charged kaons (pions) are 
needed. Since not every hadron emitted can possibly be 
detected with 100% efficiency, and since the experimen¬ 
tal techniques that identify ATs, and tt^ are very 

different, we have identified robust observables that are 
essentially independent of all these uncertainties. In par¬ 
ticular, we propose that the dynamical isospin observable 
(21) can be parameterized as in eqs. (28) and (29). DCC 
effects can appear as changes in the magnitude of the 
dynamical isospin observable as centrality is varied. We 
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emphasize that similar consequence may follow from any 
mechanism that produces many small domains that de¬ 
cay to pions and kaons, such as the Polyakov Loop Con¬ 
densate H^]. We anxiously await what RHIC will have 
to say! 
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